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Abstract 

Analytical expressions for the orbital precessions affecting the rel- 
ative motion of the components of a local binary system induced by 
Lorentz-violating Preferred Frame Effects (PFE) are explicitly com- 
puted in terms of the PPN parameters ai , a2 . No a-priori assump- 
tions about the spatial orientation of the preferred velocity w are 
made. The limit of small eccentricities is not used in the calcula- 
tion. Constraints on ai,a2 are inferred from latest determinations 
of Solar System planetary precessions. Other bounds existing in the 
literature are critically reviewed, with particular emphasis on the con- 
straint \a2\ ~ 10~^ based on an interpretation of the current close 
alignment of the Sun' s equator with the invariable plane of the So- 
lar System in terms of the action of a a2-induced torque throughout 
the entire Solar System's existence. Taken individually, the supple- 
mentary precessions Arij of Earth and Mercury, recently determined 
with the INPOPlOa ephemerides, yield ai = (0.8 ± 4) x 10"^ and 
a2 = (4 ± 6) X 10~^, respectively. A linear combination of the sup- 
plementary perihelion precessions of all the inner planets of the Solar 
System, able to remove the a-priori bias of unmodelled/mismodelled 
standard effects such as the general relativistic Lense-Thirring preces- 
sions and the classical rates due to the Sun's oblateness J2, allows to 
infer ai = (-1 ± 6) x 10-^ a2 = (-0.9 ± 3.5) x 10-^. Such bounds 
should be improved in the near future after processing the data that 
are being collected by the MESSENGER spacecraft, currently orbiting 
Mercury. Further improvements may come in the mid-future from the 
approved BepiColombo mission to Mercury. 



1 Introduction 



The invariance of the currently accepted laws of physics under Lorentz trans- 
formations of the spacetime degrees of freedom [1] is one of the most far- 
reaching ingredients of our vision of the physical reality. To date, it is rou- 
tinely corroborated with the greatest accuracy in the high-energy realms of 
particle physics [2] . The current level of confidence we have on its validity in 
gravitational physics [3] is relatively more modest because of the challenges 
in accurate experiments involving the gravitational interaction [Ij. 

Modified models of the gravitational interaction encompassing violations 
of the Lorentz symmetry arise in several theoretical scenarios [5H7| such as, 
e.g., vector-metric theories Tensor- Vector-Scalar (TeVeS) theories [9], 
Einstein- jEther theories dDHH], MOND pW8] . Hofava-Lifshitz-type theo- 
ries |19fl23j . supersymmetric field theories |24p25]. Standard Model Exten- 
sions (SME) p6H30] . due to the existence of a putative Preferred Frame 
(PF) which might be singled out by, e.g., the Universe matter distribution. 

In the framework of the Parameterized Post-Newtonian (PPN) formal- 
ism, PFE are phenomenologically taken into account by the PPN parameters 
ai,a2,a3 [8]. Since is also related to possible violations of the matter- 
energy conservation and is very accurately constrained down to a w lO"^'^ 
level from pulsar acceleration data [31] . it will be neglected in the follow- 
ing. For an overall overview on the present-day constraints on all the PPN 
parameters, see [4j. As far as ai,a2 are concerned, see also the discussion 
in Section m at present, they are constrained at « 10"'* {ai),^ 10"'' (02) 
level by a variety of techniques. 

In this paper, we intend to use latest data from planetary motions in our 
Solar System to constrain both ai and a2 to a ~ 10~^ — 10~^ level. 

The plan of the work is as follows. In Section [2] we work out the ai,a2 
Hamiltonians. They are the basis for perturbatively calculating the aver- 
aged orbital precessions of the relative motion of the binary's components in 
Section [3] by using the Lagrange planetary equations. The resulting analyt- 
ical formulas are exact in the sense that no a-priori simplifying assumptions 
about the PF velocity and the orbital geometry are made. In Section HI 
after critically discussing the constraints on qi, 02 existing in literature, we 
use latest determinations of the orbital motions of some planets of the Solar 
System to infer accurate bounds on both ai and 02- Section [5] summarizes 
our findings. 
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2 The perturbing Hamiltonians 



Let us consider a local binary system of total mass M = mA + rriB- Let ^'cM 
be the position vector of the system's center of mass (CM) with respect to 
the origin of the PF. Thus, we can write 

rl = r'p^ + rl^ = -'^r + r%^, (1) 



where the index ' refers to the CM frame; the position vector of m-Q with 
respect to is r = — and f = r/r is its unit vector. It follows 

vl = -'^v + w, (3) 



vf, = V + w, (4) 

where w is the PF velocity of the CM. 

At post-Newtonian level, the PFE two-body reducec|3 Lagrangian con- 
sists of the sum of the following terms |32H35j 



with 



^"1 = -^^K-^b), (5) 



ai 

~ . a2GM , Q / 

= (""A • r) [v^B • r) , (7) 

where G is the Newtonian constant of gravitation, and c is the speed of light 
in vacuum. For a comparison with the bounds on 02 inferred by Nordtvedt 
|32^l33j. we notice that he adopted a nonstandard normalization: a^"^'^ = 
02/2. 

Since we are interested in the relative orbital motion of the binary's com- 
ponents, we have to express eq. ([5])-eq. ([7]) in terms of relative quantities. 



^It is the ratio of the Lagrangian to the reduced mass fi = mArriB/M of the binary 
system. 
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From eq. ([Sj-eq. ^ it turns out, 



n n 2 / N 



rriAmB 2 



(8) 



« • r) (v% ■ r) 



Am 



M2 



aiGM 
2c^r 



Am 



(i? • w) 



mArriB 2 



(9) 



where Vr = v • f,Wr = w • f, Am = rriA — n-B- Thus, eq. ©-eq. become 



(10) 



with 



012 



a2GM 
2c2r 



Am 



■WrVr 



mAme 2 



M2 



fill 



(12) 



By recalling that a general velocity-dependent perturbative Lagrangian i2pert 
yields a corresponding perturbative Hamiltonian 



n 



pert 



-c 



1 /dC 



pert 



pert 



dv 



(13) 



the following PFE reduced Hamiltonians can be obtained, to order O {G/c^^ , 
from eq. ([lOj-eq. ^ 



n 



ai 



aiGM 



r, Am , , rriAm-a, 9 
w + ^r-r {v -w) TTTT^u^ 



M 



M2 



(14) 
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0-2 



with 
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ai 



a2GM 



OL2 



2 Am 



mArriB 2 



2c^r V M M2 

Perturbing accelerations can be obtained from eq. (jl4p -eq. (jl6p as 



(15) 
(16) 



by noting that 



dv 



aiGM 
2c^r 




2 



M2 
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) 



(18) 



dC. 
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02 GM 
2c^r 



2 



M2 



) 



r. 



(19) 



3 Averaged orbital perturbations 

An effective method to analyticahy calculate the orbital effects induced by 
a generic small correction ^pert to the Newtonian potential = —GM/r 
consists of evaluating Upcvt over the unperturbed Keplerian ellipse, assumed 
as reference trajectory, and averaging it over one orbital period Pb of the 
test particle. Then, the Lagrange perturbation equations |38ll39j allow to 
straightforwardly calculate to first order the long-term rates of change of 
the orbital elements by means of partial derivatives of {T-Lpert) with respect 
to them. In principle, it would be possible to adopt a different reference 
trajectory as unperturbed orbit including also general relativity at 1 PN 
level and use the so-called post Newtonian Lagrange planetary equations 
[JOlHT]. Nonetheless, in the specific case of eq. p^ -eq. (fT6]) . in addition 
to the first order precessions of order O (c~^) coming from the use of the 
Newtonian reference trajectory to be computed below, other "mixed" IPN- 
PF precessions of higher order would arise specifying the influence of PFE on 
the IPN post-Keplerian orbital motion assumed as unperturbed. From the 
point of view of constraining ai,a2 from observations, they are practically 
negligible since their magnitude is quite smaller than the first order terms. 

Before performing the calculation, a consideration about the orbital el- 
ements employed in it is in order. Strictly speaking, since the perturbing 
Hamiltonians of eq. (|14p -eq. ()16p depend also upon the velocity v, insert- 
ing them into the usual Lagrange variation equations yields nonosculating 
orbital elements called contact elements |42] . which are not tangent to the 
perturbed trajectory [36||43j . If certain general conditions on the periodicity 
of motion are satisfied [44j . the averaging procedure reduces the difference 
between the contact and the osculating elements in such a way that residual 
differences occur only to higher order. See |45] for a preliminary discussion 
of this issue in the present context. As it will be discussed in Section HI it 
is just what occurs in the present case. 

The semimajor axis a, which is dimensionally a length characterizing the 
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size of the orbit, varies according to the Lagrange equation [3^ 



da \ _ 2 dCHpevt) ^^^^ 



dt I n\^a dM. 

In it, rib = -s/ GMj ar" is the Keplerian mean motion in terms of which the 
orbital period is expressed as Pb = 27r/nb; its dimensions are those of a 
frequency. The angle M = n\^(t — to) is the mean anomaly which is one 
of the fast variables reckoning the instantaneous position along the orbital 
path; to is the time of passage at pericenter. 

The Lagrange equation for the eccentricity e, which is a dimensionless 
numerical parameter determining the shape of the ellipse in such a way thaJl 
< e < 1, is 



de\_ 1 
dt I riTaO^ 



e 



diTipert) din 



pert/ 



(21) 



In it, the argument of pericenter u is an angle in the orbital plane counted 
from the line of the nodes (see below) to the point of closest approach. 
Although not yet explicitly determined for Solar System's planets [36], the 
rate of change of e is actually measured in some binary pulsar systems 
[3511171118]. 

The inclination / of the orbital plane to the reference {x, y} plan^ Its 
rate of change has not yet been explicitly determined by astronomers for 
planets of the Solar System [46j . On the other hand, its precession is one 
of the observable quantities in binary pulsar studies through its connec- 
tion with the projected semimajor axis Xp = OpSinJ/c, where ap is the 
barycentric semimajor axis of the pulsar. Indeed, its rate of change Xp is 
often accessible to observation in several binaries 1351147^ 148]. The Lagrange 
equation for / is [38j 



dl\ 1 



dt / n\,a?'\J\ — sin/ 



d{T-Lpen) jd{l-LpcTt) 



(22) 



In it, the longitude of the ascending node J7 is an angle in the {x, y} reference 
plane counted from the x reference direction to the line of the nodes which, 
in turn, is the intersection of the orbital plane with the {x, y} plane itself. 



A circular orbit corresponds to e = 0. 



^It is the plane of the sky in the case of a celestial binary system such as, e.g., a binary 
pulsar, while it is customarily the Earth's mean equator at the epoch J2000.0 if the planets 
of our Sun are considered. 
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The node is one of the orbital parameters used by astronomers @6] 
to constrain putative non-standard dynamical effects in the Solar System. 
Sometimes, it is measurable also in specific binary pulsars |49H51j . The 
long-term precession of fl is computed from the Lagrange perturbation equa- 
tion [SB] 



The longitude of pericenter w = J7 + cj, which is a "broken" angle, is 
another parameter usually adopted in both Solar System and binary pulsars 
studies to put limits on putative exotic forces [l6]. Its Lagrange perturbation 
equation is [38] 




In principle, also the mean anomaly M should be considered, but it will 
not be treated here since its precession is neither determined in Solar System 
analyses of planetary motions nor in binary pulsars. Indeed, its accuracy is 
necessarily limited by the uncertainty in the gravitational parameter GM 
entering nb, making the use of M usually less competitive than the other 
orbital elements. 

Other analytical calculations of the orbital effects of non- vanishing ai , 02 , 
made with different computational strategies and based on various levels of 
approximation, can be found in, e.g., t34 ^ l35 t H5 t l52 y 53) . 

3.1 The tti precessions 

The orbital precessions induced by eq. ()14p can be conveniently computed 
by considering each of the three terms constituting eq. (|14p separately. 
The following analytical expressions are exact in the sense that no a priori 
assumptions on the spatial orientation of w were assumed. Moreover, the 
limits of small eccentricity and inclination were not assumed as well. 

3.1.1 The term 

The long-term precessions caused by the term proportional to the square 
of the two-body relative velocity in eq. (jl4p are all zero, apart from the 
longitude of the pericenter whose rate of change is 




(23) 





(24) 



dw 
~dt 




(25) 



7 



3.1.2 The term 

All the precessions due to the term proportional to w'^ in eq. (fT^ vanish. 

3.1.3 The {v • w) term 

The mixed term of eq. (jl4p . proportional to Am, yields non- vanishing long- 
term variations for all the orbital elements, apart from a They are 

/ de\ aiwGAm ( — 1 + + a/I — 



\di/^ 2c2a2e2 {cosa;(u)^cosr2 + u;ysin^2) + 



+ sincj [wz sin / + cos / {wy cos i7 — Wx sin /?)]} , (26) 

aiwGAm ( — 1 + + \/l — ) sin a; 
Tt) = 2cVe(l-e^) 

— sin / {wy cos Q — Wx sin J?)] , (27) 

\ aiiuGAm ( — 1 + + Vl — ) esc I cos a; 
-dt) = 2c^a^e(l-e^) ^"'^ 

— sin / {wy cos n — Wx sin J?)] , (28) 



aiwGAm I / — l + Vl — 



[ri;^ sin /cos 0;+ 



+ cos / cos oj {wy cos Q — Wx sin J?) 
— sin a; {wx cos fl + Wy sin /?)] — 



-l + e^ + Vr^\, (I\ J 

tan - cos UJ [w-r cos 7 - 



l-e2 i V2 
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— sin I {wy COS n — Wx sin n)]} . (29) 

It turns out that the leading term in eq. ()26p is of order O (e^) in the 
eccentricity, while the next-to- leading order one is of order O (e^). As far as 
/ and n are concerned, their leading order terms are of order O (e), while 
the first non-vanishing terms of higher order in e are of order O (e^) . The 
pericenter precession of eq. (j29p has a term of order O (e~^) , while the next 
ones are of order O (e) and O (e'^) , respectively. 

3.2 The a2 precessions 

The orbital precessions due to Q2 consist of the same as of Section 13.11 
rescaled by —02/0:1, and by those coming from eq. (fT6]l . Also in this case, 
analytical expressions exact in w, e, / are obtained. 

3.2.1 The term 

All the long-term precessions induced by the term proportional to in eq. 
(jl6p vanish, apart from zu. Its rate is 

dw \ a2GmpjnBn\, ^^^^ 



dt I 2c^a{\-e^)M' 
3.2.2 The wl term 

The non- vanishing long-term rates of change of the orbital elements due to 
wl in eq. (fTHl) are as follows. 



X 



\ a2'«;2nbVl - (-2 + + 2Vl - e 
~dtl~ 8c2e3 

X \^—%Wz sin / cos 2bj (w^ cos H + Wy sin Q) + 

+ 4 cos / cos 2a; \—2wx'Wy cos 2Q + [wl. — w^) sin 2ff\ + 

+ sin 2oj [[wl - wl) (3 + cos 21) cos 2/?+ 

+ 2 sin / [w^ + Wy — 2w^) — Awz sin 21 [Wy cos Q — Wx sin Q) + 
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+ 2wrcWy (3 + COS 21) sm2n] } , 



dl\ a2W-'^"'b 
di/ ~ 2^^VT^ 



[wz COS / — sin / (wy cos Q — Wx sin Q)] x 



X |— 2\/ 1 — e^Wx cos Q s\r? uj+ 

+ 2\/l — e^Wy cos / cos fi sin 2a;— 

— 2\/l — e^Wy sin sin^ u+ 

+ 2\/ 1 — cos^ a; (i&a; cos Q + Wy sin 



— 2 cos 2(jj {wx cos n + Wy sin i?) — 



2e sin a; {wx cos + sin /?) + 



+ 2\/r^-^ sin 2a; [iD;^ sin I — Wx cos J sin fi] + 

+ (—2 + e^) sin 2a; [wz sin 7 + cos I {wy cos fi — Wx sin i?)] } , 



dt I 2c2e2Vr^ 



\wz COS J — sin 7 (u;^^ cos fi — Wx sin x 



X ^Wz sin/ + 



+ f -2 + + 2^1 - e^) cos 2a) 



+ 



+ cos/ 



+ (|-2 + + 2V1 - 



cos 2a; 



{wy COS fi — Wx sin 



— ^— 2 + + 2-\/ 1 — e^^ sin 2a; cos fi + Wy sin i?) | , 



10 



dt / 8c2e4 



{(-2 + e2 + 2\/l-e2) 



[( 



X \{wl - Wy) (3 + COS 21) cos 2 J7 cos 2uj+ 



+ 2{wl + wl- 2wl) sin^ I cos 2uj- 



sin 2/ cos il — Wx sin /?) cos 2a;+ 



+ 2wxWy (3 + cos 2/) sin 2Q cos 2a;+ 



+ Stf)^ sin / {wx cos i7 + sin Q) sin 20;+ 



+ St^ajWu cos /cos 2i7 sin2a;- 



— 4 (w^ — cos / sin 2i? sin 2a;] 



4e2tan (^) . ^ , „ ^ . ^x-, 

— [w^ cos / — sin I [Wy cos U — Wx sin J7)J x 



^+ (-2 + e2 + 2\/l-e2) cos2w) sin/+ 



+ cos I ie^ + i —2 + + 2\/ 1 — ) cos 2a; ) (wj^ cos Q — Wx sin J?) 



(^-2 + + 2\/l - (w^ cos i? + sin i?) sin 2a; | . (35) 



In the Umit e — )• 0, the leading terms in eq. (j3T]) -eq. (p5|) are of or- 
der O (e'^). The next non- vanishing terms are of order O (e^) for ah the 
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elements, with the exception of e whose next-to- leading term is of order 
0(e). 

3.2.3 The VrWr term 

The mixed term of eq. (jl6p . proportional to Am, causes non- vanishing long- 
term changes for all the orbital elements, apart from a. They are formally 
identical to those of Section r3.1.3l with ai replaced by 02- Thus, Am does 
not contribute to the overall orbital perturbations due to 02 because the 
precessions coming from the Am term in eq. p6p are canceled by those 
coming from the Am term in the first piece of eq. (jlSp . 

The total 02— induced precessions are made up of the {w/c)^ terms of 
Section [3.2.21 proportional to a~^/^, and of the terms containing the product 
of the masses which are completely negligible in the Solar System. In the 
case of a tight binary pulsar with, say, ~ ni-B ~ M/2, a ~ 9.0 x 10^ km, 
the precessions of Section [3.1.11 rescaled by — a2/aii and Section [3.2.11 are 
proportional to (vorb/c)^, with Vorh ~ 600 km s""*^; for typical values of w 
(see Section HI, they are dominant with respect to the precessions of Section 

4 Confrontation with the observations 

The analytical results of Section 13.11 and Section 13.21 were successfully con- 
firmed by numerically integrating the equations of motion of a test particle 
including the PF accelerations coming from eq. ()17p with eq. ()14p -eq. ()16p 
and eq. ()18p -eq. ()19p over a given time interval, and plotting the temporal 
evolution of the numerically computed orbital elements with the standard 
formulas of the osculating Keplerian elements. In Figure [T] and Figure [2] 
the plots for Mercury over one century, obtained for qi = 02 = 1 just for 
illustrative purposes and with eq. ()54p for are displayed. They yield 
the same centennial rates obtainable from the analytical formulas of Section 
13.11 and Section 13.21 for the contact elements, computed for the same values 
of ai,a2,w used in the numerical integration. This shows that, to the ap- 
proximation level used in the calculation, the contact elements are equal to 
the osculating elements. 

Thus, confident in our results of Section [STTl Section 13.21 we are going to 
apply them to effectively constrain both ai and 02 from latest determina- 
tions of planetary orbital motions |46] . In Section [4.11 we critically review 
the bounds existing in literature. 
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Figure 1: Numerically integrated temporal evolutions of the total PF 
ai— induced shifts (in arcsec) of the eccentricity e, the inclination /, the 
longitude of the ascending node H and the longitude of perihelion vj of 
Mercury over one century. They were obtained by numerically integrat- 
ing the planet's equations of motion in cartesian coordinates with and 
without the perturbative accelerations arising from eq. p4|) by using the 
same initial conditions retrieved from the NASA JPL HORIZONS Sys- 
tem (http://ssd.jpl.nasa.gov/7horizons). For illustrative purposes the value 
ai = 1 was chosen. The figures in eq. (jSl]) . taken from [5l], were adopted 
for w. The resulting centennial rates are in agreement with the sum of those 
computed from the analytical formulas in Section 13.11 
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PF a2 shift of Mercury's eccentricily (iTt - 1) PF o-t shift of Mercury's inclination (o'2 - 1) 




Figure 2: Numerically integrated temporal evolutions of the total PF 
0(2— induced shifts (in arcsec) of the eccentricity e, the inclination /, the 
longitude of the ascending node H and the longitude of perihelion uj of 
Mercury over one century. They were obtained by numerically integrating 
the planet's equations of motion in cartesian coordinates with and with- 
out the perturbative accelerations arising from eq. P^ -eq. by using 
the same initial conditions retrieved from the NASA JPL HORIZONS Sys- 
tem (http://ssd.jpl.nasa.gov/7horizons). For illustrative purposes the value 
a2 = I was chosen. The figures in eq. (jSl]) . taken from [5l], were adopted 
for w. The resulting centennial rates are in agreement with the sum of those 
computed from the analytical formulas in Section 13.21 
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4.1 Discussion of the constraints existing in literature 

As a general rule pertaining the bounds on the strong- field equivalent ai, 02 
with compact object, Shao and Wex [35] interestingly remark that their 
general validity may be questioned by a potential compactness-dependence 
(or mass-dependence) of ai, 02 because of certain peculiar phenomena, such 
as spontaneous scalarization [55j, which may take place. Thus, Shao and 
Wex [35] warn that it is always recommendable to specify the binary system 
used to infer given constraints on ai,a2- Another limitation of the strong- 
field tests is their probabilistic nature since the orbital configurations of the 
binaries used are not always completely known: indeed, the longitudes of 
the ascending nodes J7 are often unknown, and also the inclinations / may 
be sometimes determined modulo the ambiguity of / — )• 180° — /. 

4.1.1 ai 

Earlier constraints on ai were obtained from Solar System planetary motions 
[52|[55] and Lunar Laser Ranging ^57j ; the most accurate bounds [56n57j were 
at 2 X 10-^ level. 

Strong-field tests were reported by using compact objects such as neutron 
stars and white dwarfs [331 EH]. They were able to constrain di, which, 
in principle, may differ from the weak-field value ai because of possible 
contributions from strong-field effects, down to 5.0 — 1.2 x 10~^. Latest 
bounds on qi were obtained by Shao and Wex [35] from the neutron star- 
white dwarf binary PSR J1738+0333 [59] yielding |di| = -0.4l|{ x 10"^. 

4.1.2 a2 

Concerning Q2, Nordtvedt |33j inferred an upper bound of the order of 10^^ 
from the close alignment of the Sun's angular momentum with the total 
angular momentum of the Solar System. Nonetheless, as remarked by Shao 
and Wex [35], such a result strongly depends on the assumption that the two 
angular momenta were aligned just after the formation of the Solar System 
about 5 Gyr ago. If compared to the well known planetary orbital motions 
and the dynamical forces determining them, the hypothesis on which the 
analysis by Nordtvedt [33] is based is relatively more speculative and less 
testable. Another potential issue of the Nordtvedt reasoning [33j may be 
the following one. 

The orientation of the Solar System's invariable plane |60) is fixed by 
the values of the celestial coordinates of its north pole: the right ascension 
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a (RA) and the declination 5 (DEC). At epoch J2000.0, they are [6l] 

ainv = 273.85°, (5inv = 66.99°; (36) 

thus, its normal unit vector L is 

Lx = cos Jinv cos ainv = 0.03, (37) 

Ly = cos 5inv sin ainv = —0.39, (38) 

= sin (5inv = 0.92. (39) 

The north pole of rotation of the Sun at J2000.0 is characterized by [61] 

a© = 286.13°, 5© = 63.87°, (40) 
so that the Sun's spin axis is 

5® = 0.12, (41) 



5® = -0.43, (42) 



S® = 0.89. (43) 

Thus, the angle between the Sun's spin axis S and the Solar System's total 
angular momentum L is 

Ols = 5.97°. (44) 

While the orientation of the invariable plane is nowadays quite accurate, 
being at « 0.1 — 0.3 mas cty~^ level [60], it is not so for the Sun's spin 
axis. Indeed, its determination is usually made in terms of the Carrington 
elements [621163] i, which is the angle from the Sun's equator to the ecliptic, 
and the longitude of the node of the Sun's equator with respect to the 
Vernal equinox Y along the ecliptic, by means of 

s® = sinisinri, (45) 
Sy=— sin i cos fi, (46) 
s® = cosi. (47) 
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Beck and Giles |64) recently measured them from time-distance helioseis- 
mology analysis of Dopplergrams from the Michelson Doppler Imager (MDI) 
instrument on board the SOlar and Heliospheric Observatory (SOHO) space- 
craft; the data span was At = 5 yr, from May 1996 through July 2001. The 
resulting values ar^ 

i = 7.155° ±0.002°, (48) 



n = 73.5° ± 1°. (49) 

Maybe a better accuracy could be reached in a near future with a dynam- 
ical measurement from planetary orbital motions by exploiting the general 
relativistic Lense-Thirring effect p5j. The figures of eq. (|48p -eq. (|49|) 
imply an ability to observationally constrain putative rates of change of 
Ols {i,^ , Qinv; '^inv) with a necessarily limited accuracy; it can approximately 
be evaluated as 



(jg^^ — — ^ = 7.73 X 10 arcsec yr. (50) 

It is an important limiting factor when hypothesized secular rates 6ls 
used to constrain parameters entering the models proposed to explain Oi^g. 

It is just the far as 02 is concerned. Indeed, Nordtvedt [33] 

noticed that one of the dynamical consequences of the first term in eq. (jl6p 
is a torque causing a precession of the Sun's spin axis about w at a, rate 
proportional to 02- Competing classical torques due to the planets of 
the Solar System induce an extremely low overall precession of the Sun's spin 
axis characterized by a rate [33] Hdass ~ 10^^*^ yr. Under certain simplifying 
assumptions, Nordtvedt [33] obtained 



\ / \ ■ ( 'class^ \ 



sm I I = I „ I sm 



- 1 ■ (51) 



By assuming that the Sun's equator and the invariable plane were aligned at 
the birth of the Solar System 5 Gyr ago, and that the present-day value of eq. 
(dD) is due to the steady action of the aforementioned torques throughout 
the life of the Solar System in such a way that HciassAT/2 <C 1, AT = 5 Gyr, 
eq. ([5T]) reduces to 

Ols ~ H„,Ar. (52) 



^ After rotating eq. (|45p -ea. (I47p . calculated with eq. (|48[) -ea. (|49l) . from the ecliptic 



to the Earth's equator, eq. (|41l) -ea. (|43p are obtained. 
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In other words, the PF 02— induced spin precession would cause a secular 
rate of change Ols whose magnitude is approximately equal to H^j. By 
posing AT = 5 Gyr, Nordtvedt [33j inferred an upper limit of a2 as little as 
w 10~^ from a hypothesized precession 

5 97° 

e^s = — — = 4 X 10"^ arcsec yr"\ (53) 
5 Gyr 

In fact, apart from the more or less speculative assumptions about a dy- 
namical evolution spanning 5 Gyr, the observational uncertainty in actually 
measuring such a kind of time derivative of O^s should be considered in de- 
riving realistic constraints on 02- As we have seen, the errors in determining 
the Carrington elements of the Sun's spin axis yield an uncertainty in ob- 
servationally constraining Ols which may be about 8 orders of magnitude 
larger than eq. (153]) . 

For previous tests involving also other phenomena such as, e.g., Earth 
tides, see [3|[66]. 

Moving to the strong-field regime, Shao and Wex [35| reported |q!2| < 
1.8 X 10~^ from the rate of change of the projected semimajor axis of the two 
pulsar-white dwarf binaries PSR J1012+5307 [67] and PSR J1738+0333 [59]. 
Earlier limits (—0.3 < a2 < 0.2) were inferred by Wex and Kramer [53| 
by using the double pulsar PSR J0737-3039A/B [68] made of two neutron 
stars. As remarked by Shao and Wex [35] themselves, a potential drawback 
of their latest constraint on 0:2 is that it was obtained by combining data 
from two systems whose neutron stars have different masses. Moreover, 
a straightforward comparison with the previous test by Wex and Kramer 
[53] is, in principle, problematic because the latter probed the interaction 
between two strongly self-gravitating objects, while the systems used by 
Shao and Wex 1351 contain white dwarfs. 



4.2 Bounds from recent Solar System planetary data 

According to Table [H the present-day accuracy in constraining possible 
extra-precessions of the nodes and perihelia with respect to standard New- 
tonian/Einsteinian effects has reached the ~ 1 — 0.5 milliarcseconds per 
century (mas cty~^) level for the inner planets of the Solar System and for 
Saturn. This fact, in conjunction with the theoretical predictions of Section 
[31 allows to infer tight constraints on both ai and a2- 

As far as w is concerned, a natural choice for the preferred frame, com- 
mon to the literature on preferred- frame effects [33H35| H5 | [56 l [66]. is the 
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Table 1: Supplementary precessions Ai?, Aro of the longitudes of the node 
and of the perihelion for some planets of the Solar System estimated by 
Fienga et al. |46) with the INPOPlOa ephemerides. The reference {x,y} 
plane is the mean Earth's equator at J2000.0. The units are milliarcseconds 
per century (mas cty~^). 





A/2 (mas cty ^) 


Azu (mas cty ^) 


Mercury 


1.4 ±1.8 


0.4 ±0.6 


Venus 


0.2 ± 1.5 


0.2 ± 1.5 


Earth 


0.0 ±0.9 


-0.2 ±0.9 


Mars 


-0.05 ±0.13 


-0.04 ±0.15 


Saturn 


-0.1 ±0.4 


0.15 ±0.65 



Cosmic Microwave Background (CMB). With such a choice one is quite plau- 
sibly assuming that the preferred frame is determined by the global matter 
distribution of the Universe, with the extra-components of the gravitational 
interaction ranging over scales at least comparable to the Hubble radius. 
Latest results from the Wilkinson Microwave Anisotropy Probe (WMAP) 
yield a peculiar velocity of the Solar System Barycenter (SSB) of [Si] 

w = 369.0 ± 0.9 km s-\ / = 263.99° ± 0.14°, 6 = 48.26° ± 0.03°, (54) 

where I and b are the Galactic longitude and latitude, respectively. 

By comparing the figures in Table [1] with the theoretical predictions of 
Section [3l computed with eq. ([5l|) . it is possible to infer upper bounds on 
ai and 02- It turns out that, in general, the nodes yield weaker constraints 
than the perihelia, especially as far as ai is concerned. Thus, in Table[2]only 
the bounds from the perihelia are displayed. The most stringent constraint 
on ai comes from the perihelion of Earth, amounting to |ai| ^ 10~^. The 
perihelion of Mercury yields the tightest constraint for 02, which is of the 
order of \a2\ ^ 10~^ as well. In principle, the bounds of Table [2] may 
be biased by certain standard dynamical effects which were not modeled 
at all by Fienga et al. [l6], such as the general relativistic Lense-Thirring 
effect [69] , thus impacting the supplementary rates Aw of Table [H The 
gravitomagnetic precession should not be neglected since its expected value 
(—2.0 mas cty""*^) is larger than the current uncertainty in the Mercury's 
perihelion extra-rate (0.6 mas cty~^); see [65] for a recent discussion of this 
aspect. Another standard effect which may a-priori impact the bounds in 
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Table 2: Constraints on ai and 02 obtained from a straightforward compar- 
ison of the figures of Table [T] for the supplementary rates of the planetary 
perihelia with the theoretical predictions of Section O calculated with eq. 
(j54p . The tightest bounds comes from the perihelion of Earth (ai) and of 
Mercury (02)- 







02 


Mercury 


(-3 ±5) X 10-6 


(4 ± 6) X 10-6 


Venus 


(-0.1 ± 1.1) X 10"^ 


(-0.7 ±5.7) X 10"^ 


Earth 


(0.8 ±4) X 10"*^ 


(-0.8 ±3.7) X 10-5 


Mars 


(-0.8 ±2.9) X 10"^ 


(0.4 ± 1.5) X 10-5 


Saturn 


(-1.94 ±8.41) X 10"^ 


(1.95 ±8.5) X 10-"^ 



Table [2] through its lurking in Table [His the Newtonian perihelion precession 
due to the Sun's quadrupole moment J2. Although its dynamical effect was 
fully modeled by Fienga et al. [46], a lingering uncertainty of [70] ~ 10% still 
affects J2 in such a way that a resulting mismodelled perihelion precession 
certainly contributes to the figures in Table [H The availability of more than 
one supplementary perihelion rate is a great advantage since it allows to set 
up a suitable linear combination involving the perihelia of the four inner 
planets able to separate, by construction, ai, 02 from J2 and the Lense- 
Thirring effects. From the following linear system of four equations in the 
four unknowns ai, 02, J2, fJ-LT 

Auj^ = aiw^a^ + a2W^,a2 + '^'^'^^2 t^^'T^^^T^ j = Mercury . . . Mars, (55) 

where the coefficients w, are the analytical expressions of the pericenter 
precessions caused by the effects considerecil, it is possible to obtain 

ai = (-1±6) X 10-6, (56) 

02 = (-0.9 ±3.5) X 10-5. (57) 

By construction, such values are not a-priori biased by any mismodeling in 
J2 and by the Lense-Thirring effect. 

The constraints of Table [2] and of eq. (f56l) -eq. ([57|) will likely be fur- 
ther improved in the next few years. Indeed, our knowledge of the orbit of 

^General expressions for the J2 and Lense-Thirring pericenter precessions can be found 
in [71]. The coefficient /ilt is one in general relativity. 
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Mercury will be greatly enhanced when the entire data set of the MESSEN- 
GER spacecraft |72j . which was inserted in orbit around Mercury in March 
2011 for a year-long science phase extended until March 2013, will be pro- 
cessed. Another planned mission to Mercury is BepiColombo [73j: it will be 
launched in 2015 and should reach its target in 2022. One of its experiment 
is MORE [71]; it should further improve the determination of the orbit of 
Mercury [75] . 

5 Summary and conclusions 

The long-term, i.e. averaged over one orbital revolution, precessions of the 
orbital elements of the relative motion of the components of a binary system 
due to non-zero values of the PPN parameters ai and 02 inducing preferred- 
frame effects were analytically worked out to first-order with the perturba- 
tive Lagrange planetary equations. We did not recur to a-priori simplifying 
assumptions on the preferred-frame velocity w and on the binary's orbital 
geometry. The analytical expressions of the precessions were successfully 
tested with a numerical integration which yielded the same results. It was 
also numerically checked that, to the order of approximation used, the con- 
tact orbital elements used in the analytical calculation are actually tangent, 
i.e. osculating, to the perturbed trajectory. 

We critically discussed the existing bounds on 01,02- We remarked that 
the O (10^^) constraint on 02 inferred by Nordtvedt from a hypothesized 
secular precession of the angle between the Sun's equator and the invariable 
plane of the Solar System throughout its entire existence should be consid- 
ered as optimistic. Indeed, it did not take into account the measurement 
limits in actually constraining such a precession from observations in view 
of the lingering uncertainties in the Carrington elements i and, especially, 17 
which determine the orientation of the Sun's spin axis. 

We used the supplementary node and perihelion secular precessions AQ, Aw 
of some planets of the Solar System, recently determined by astronomers 
from an analysis of a centennial data record in building the INPOPlOa 
ephemerides, to effectively constrain both oi and 02- Indeed, such supple- 
mentary precessions account for any unmodelled/mismodelled dynamical 
effects with respect to the standard Einstein-Newton gravity. After having 
noticed that the perihelia give tighter constraints than the nodes, we lin- 
early combined the supplementary perihelion precessions Aw of Mercury, 
Venus, Earth, Mars to simultaneously solve for 01,02, and for the Sun's 
Lense-Thirring effect and oblateness J2 as well. Indeed, they were unmod- 
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elled/mismodelled in the INPOPlOa-based analysis and represent competing 
effects for tlie ai, 02— induced precessions we are interested in. We obtained 
ai = (-1 ± 6) X 10-^a2 = (-0.9 ± 3.5) x 10"^ Our results, which retain 
a general validity for all the weak- field scenarios, are based on well tested 
dynamical effects pertaining known Solar System's objects, and do not rely 
upon speculative assumptions about hypothesized phenomena for which no 
data actually exist. 

Future improvements in the orbit determination of Mercury are 
expected in the near and mid future from MESSENGER and BepiColombo 
missions; they should allow to get tighter constraints on 01,02. 
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